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DESCRIPTION OF THE INVENTION 

Related Applications 

[001] This application claims priority to U.S. Provisional Application Serial 
No. 60/444,281 , filed on January 31 , 2003, which is incorporated herein by 
reference. 

Field of the Invention 

[002] The invention relates generally to methods for reducing the time 
and/or cost associated with extraction of hydrocarbons from underground reservoirs. 
More specifically, the present invention relates to systems and methods for 
automating the selection of targets for well placement, using two-dimensional 
matrixes that represent a three-dimensional model of the reservoir. 

Background of the Invention 

[003] Historically, to amve at a set of targets for the placement of wells for 
the extraction of hydrocarbons, a manual process is used in which the surfaces 
representing the top, bottom surface, or intermediate surface of the reservoir is 
interpreted from seismic data to place the targets. In this process, the surface is 
overlaid with a color-coded attribute overlay, which may be comprised of attributes 
indicating the porosity, permeability, saturations, or other characteristics of the 
reservoir. Based on the color-coded attributes, the target locations are visually 
selected, the attribute color coding being the primary indicator of the target location. 

[004] In the above method, since the color coding is used as the primary 
cue for target location, it is very difficult to take into consideration the drainage radius 
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and an entire three-dimensional modei of a reservoir, because only the attributes of 
a surface are being considered. To take into consideration the entire three- 
dimensional model of the reservoir, this method would entail sequentially evaluating 
adjacent layers of the reservoir until the entire reservoir is completely evaluated, 
while attempting to remember the attributes associated with each of the previous 
layers. This a very tedious and time consuming process. 

[005] In another historically used and more accurate method than the 
above, the raw seismic data is visually analyzed for so called "bright spots." The 
bright spots are basically high or low amplitude attribute areas indicative of the 
porosity, pemneability, saturations, or other characteristic of the reservoir. In this 
process, a geocellular model of the entire reservoir is created rather than evaluating 
a selected horizontal surface of the reservoir and filtering the model to only display 
the areas or "bright spots" that have certain desirable attributes. In similar fashion to 
the previous method, in this method, the targets are also manually selected on the 
resulting geobodies. This is also a time consuming and inaccurate method of 
locating targets. 

[006] Therefore, there is a need for an automated method for locating 
targets for a hydrocarbon reservoir. A method where a user sets a series of filters 
for processing a three-dimensional geocellular model, an expected drainage radius 
for each target, and the maximum number of allowable targets. Then a computer 
implemented method traverses the geocellular model targeting areas that meet the 
filtering criteria and target spacing constraints. The method may eliminate 
overiapping targets; therefore, the method may optimize the location of each target 
by optimizing the cells that are targeted. If more target areas meet the filtering and 
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spacing criteria than the maximum allowable number, the method may only display 
targets in the most desirable areas by evaluating the cells under each target both 
horizontally and vertically and weighing each target location. 

[007] This automated process would be more desirable than the process 
discussed above of interpreting a surface of the reservoir, because the entire three- 
dimensional model of the reservoir may be taken into account, and the implemented 
selection process results in a more accurate selection of desirable target locations 
than the human selections based on visual inspection of the attributes overlaid on 
the surface. Furthermore, this automated process would be more accurate when 
compared to the method using the geocellular model, because the spacing radius 
may be optimized and the targets are weighted based on the cells under each target 
so that the placement is much more optimized than if done by human analysis. This 
automated process may result in an increase in accuracy and substantial time 
savings in the process of locating targets for well placement. 

[008] Embodiments of the present invention are directed at overcoming one 
or more of the above deficiencies described in the art. 

Summary Of The Invention 

[009] In accordance with an exemplary embodiment of the present 
invention, methods and systems are provided for automated reservoir targeting. For 
automated reservoir targeting, the systems implement methods comprising filtering a 
plurality of values in a three-dimensional model to eliminate values which are below 
a threshold to create a filtered three-dimensional model; developing a first matrix 
from the filtered three-dimensional model representing a two-dimensional model of 
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the reservoir, wherein the first matrix includes a plurality of cell center locations, cell 
areas, and the plurality of values; developing a second matrix from the first matrix 
using a distance-weighted sum of the plurality of values; and selecting target 
locations from the second matrix based on the distance-weighted sum of the plurality 
of values. 

[010] Additional objects and advantages of the invention will be set forth in 
part in the description which follows, and in part will be obvious from the description, 
or may be learned by practice of the invention. The objects and advantages of the 
invention will be realized and attained by means of the elements and combinations 
particularly pointed out in the appended claims. 

[01 1] It is to be understood that both the foregoing general description and 
the following detailed description are exemplary and explanatory only and are not 
restrictive of the invention, as claimed. 

[012] The accompanying drawings, which are incorporated in and 
constitute a part of this specification, illustrate several embodiments of the invention 
and together with the description, serve to explain the principles of the invention. 

Brief Description Of The Drawings 

[013] Figure 1 is a system environment in accordance with the principles of 
the present invention. 

[014] Figure 2 illustrates an exemplary three-dimensional model 
representing the distribution of porosity associated with an exemplary reservoir in 
accordance with the principles of the present invention. 
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[015] Figure 3 illustrates an exemplary three-dimensional model 
representing the distribution of oil saturation associated with an exemplary reservoir 
in accordance with the principles of the present invention. 

[016] Figure 4 illustrates an exemplary schematic model including a column 
of exemplary cells in accordance with the principles of the present invention. 

[017] Figure 5 illustrates a three-dimensional model representation of a 
filtered three-dimensional model in accordance with the principles of the present 
invention. 

[018] Figure 6 illustrates a sliding window applied to an exemplary column 
of cells in accordance with the principles of the present invention. 

[019] Figure 7 illustrates the location of a first target location on a three- 
dimensional model of a reservoir in accordance with the principles of present 
invention. 

[020] Figure 8 illustrates the location of a second target on a three- 
dimensional model of a reservoir in accordance with the principles of the present 
invention. 

[021] Figure 9 illustrates all the located targets on a three-dimensional 
model of a reservoir in accordance with the principles of the present invention. 

[022] Figure 10 is a flow chart illustrating the stages associated with the 
targeting method in accordance with the principles of the present invention. 

Description Of The Embodiments 

[023] Reference will now be made in detail to the exemplary embodiments 
of the invention, which are illustrated in the accompanying drawings. Wherever 
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possible, the same reference numbers will be used throughout the drawings to refer 
to the same or like parts. 

System Architecture 

[024] By way of a non-limiting example, Figure 1 illustrates a computer 
system in which the features and principles of the present invention may be 
implemented. As illustrated in the block diagram of Figure 1, a system environment 
consistent with an embodiment of the present invention may include an input module 
110, an output module 120, a computing platform 130, and a database or file system 
140. Computing platform 130 is adapted to include the necessary functionality and 
computing capabilities to implement the automated target select and platform 
generation methodology through the associated components (input module 110, 
output module 120, and database or file system 140). 

[025] In the embodiment of Figure 1 , computing platform 130 may comprise 
a PC or PDA for performing various functions and operations of the invention. 
Computing platform 130 may be implemented, for example, by a general purpose 
computer selectively activated or reconfigured by a computer program stored in the 
computer, or may be a specially constructed computing platform for carrying out the 
features and operations of the present invention. Computing platform 130 may also 
be implemented or provided with a wide variety of components or subsystems 
including, for example, one or more of the following: one or more central processing 
units, a co-processor, memory, registers, and other data processing devices and 
subsystems. Computing platform 130 also communicates or transfers dynamic 
analysis input and output to and from input module 110 and output module 120 
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through the use of direct connections or conrtmunication links, as illustrated in Figure 
1. 

[026] Alternatively, comnfiunication between computing platform 1 30 and 
modules 1 10, 120 can be achieved through the use of a network architecture (not 
shown). In the alternative embodiment (not shown), the network architecture may 
comprise, alone, or in any suitable combination, a telephone-based network (such 
as a PBX or POTS), a local area network (LAN), a wide area network (WAN), a 
dedicated intranet, and/or the Internet. Further, it may comprise any suitable 
combination of wired and/or wireless components and systems. By using dedicated 
communication links or shared network architecture, computing platform 130 may be 
located in the same location or at a geographically distant location from input module 
110 and/or output module 120. 

[027] Input module 1 10 of the system environment shown in Figure 1 may 
be implemented with a wide variety of devices to receive and/or provide the data as 
input to computing platform 130. As illustrated in Figure 1 , input module 1 10 
includes an input device 1 1 1 , a storage device 112, and/or a network 113. Input 
device 1 1 1 may include a keyboard, a mouse, a disk drive, video camera, magnetic 
card reader, or any other suitable input device for providing customer information to 
computing platform 130. Memory device may be implemented with various forms of 
memory or storage devices, such as read-only memory (ROM) devices and random 
access memory (RAM) devices. Storage device 1 12 may include a memory tape or 
disk drive for reading and providing information on a storage tape or disk as input to 
computing platform 120. Input module 110 may also include network interface 113, 
as illustrated in Figure 1 , to receive data over a network (such as a LAN, WAN, 
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intranet or the Internet) and to provide the same as input to computing platform 130. 
For example, network Interface 113 may be connected to a public or private 
database over a network for the purpose of receiving information about the 
customers from computing platform 130. 

[028] As illustrated in Figure 1 , output module 120 includes a display 
adapter 121 , a printer device adapter 122, and/or a network interface 123 for 
receiving the results provided as output from computing module 120. The output 
from computing platform 130 may be displayed or viewed through display adapter 
121 (such as a CRT or LCD) and printer device adapter 122. If needed, network 
interface 123 may also be provided to facilitate the communication of the results 
from computer platform 130 over a network (such as a LAN, WAN, intranet or the 
Intemet) to remote or distant locations for further analysis or viewing. 

Automated Target Selection 
Operational Description 

[029] In an exemplary embodiment of the present invention, computer 
automated methods utilize a three-dimensional geocellular model of the reservoir. A 
user sets a series of filters to be applied to the three-dimensional geocellular model, 
an expected drainage radius for each target, and the maximum number of allowable 
targets. Then the automated methods traverse the geocellular model targeting 
areas which meet the filtering criteria and the predetermined target spacing 
constraints. If overlapping targets are identified, the overlapping targets are 
eliminated so the methods optimize the location of each target by optimizing the 
desirable cells that are targeted. If more target areas meet the filtering and spacing 
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criteria than the maximum allowable number, the algorithm may display targets in 
the most desirable areas by evaluating the cells under each target both horizontally 
and vertically and weighing each target location. 

[030] In an initial stage of the methods of one embodiment, seismic data 
sets are used to build a three-dimensional model. Geoscientists use a modeling 
program, such as Decision Space™ Powermodel™, to interpret the seismic data 
sets both for the top and bottom of the reservoir in question to develop the three- 
dimensional model. PowerModel™ is built on gOcad technology developed by the 
gOcad Consortium and is available from Landmark, a Halliburton company, of 
Houston, Texas. The three-dimensional model provides a representation of the 
reservoir underneath the surface. The characteristics of the three-dimensional 
model, e.g., porosity, may be indicative of whether there are hydrocarbons in the 
reservoir. 

[031 ] In addition to the seismic data, exploratory drilling may also be 
utilized along with logging tools to retrieve readings that are then correlated to the 
seismic data sets. For example, measurements may be obtained for gamma ray, 
neutron porosity, or bulk density. Once the location of these measurements are 
determined, the measurements are correlated to the seismic data sets to improve 
the three-dimensional model. The resulting three-dimensional models are a set of 
irregularly spaced and perhaps irregularly shaped six sided cells with eight comer 
points and one or more attributes stored within. All of these cells may be combined 
to generate a three-dimensional model of a particular reservoir. 

[032] Figures 2 and 3, respectively, illustrate an exemplary three- 
dimensional model representing the distribution of the porosity and oil saturation 
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values associated with an exemplary reservoir. In Figures 2 and 3, the represented 
values are shown on a gray scale. The dark shade areas (202 and 302) represent 
areas containing high concentrations of the represented values (porosity or oil 
saturation) and the light shaded areas (204 and 304) represent areas containing low 
concentrations of the represented value. The areas with concentration values 
ranging between the high and low concentrations are shown in a gray scale 
proportionate to the level of concentration. 

[033] Once the three-dimensional model is generated, targets are generally 
selected based on a particular property or attribute, while making sure that certain 
other properties are within certain ranges. The targets may also be selected to 
obtain proper spacing. Target spacing is a consideration in the selection process 
because the location of the particular property of interest, e.g., porosity or oil 
saturation, may be concentrated in a specific area of the reservoir, and it is 
undesirable to locate multiple targets within a concentrated area that can possibly be 
served by one target. The target selection process will be discussed in greater detail 
below. 

[034] In the method of one exemplary embodiment of the invention, the 
three-dimensional model may be triangulated to create a three dimensional grid (X, 
Y, and Z axis) of cells (schematic model). Figure 4 illustrates an exemplary 
schematic model including a column 400 of exemplary cells (402-412). Once the 
triangulation process is complete, for each X and Y coordinate, the schematic model 
is traversed in the Z direction to filter out cell locations that do not meet one or more 
predetermined selection criteria (e.g., porosity > .2 and/or oil saturation > .3). 
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[035] Figure 5 illustrates a three-dimensional representation of a filtered 
three-dimensional model. This filtered model may correlate to the three-dimensional 
models of Figures 2 and 3. 

[036] Once the three dimensional model has been filtered, the preferred 
method may determine the most desirable set of contiguous cells that meet one or 
more of the filtering criteria. The most desirables cells may be located using a 
windowing technique in which a sliding window composed, for example, of three 
cells or three layers may be used to find the location of the most desirable three 
continuous layers. The window size may be determined by the number of layers in 
the three-dimensional model. Figure 6 provides an illustrafion of a sliding window 
being applied to cells 402-406. 

[037] If cells 402-406 are determined to be the most desirable cells in 
column 400, then the location of the center point or the midpoint 602 of these three 
layers is identified along with the value of the particular property of interest, for 
example the porosity or oil saturation value associated with the center point. The 
forgoing process distills the schematic model down to a two dimensional value 
matrix representing the most desirable layers in the columns of the schematic 
model. Every column in the model is represented by row and column indices (matrix 
entry) in the value matrix and each matrix entry may contain the value of the 
property of interest (e.g., porosity or oil saturation). 

[038] Table 1 below illustrates an exemplary value matrix for a distilled three- 
dimensional model. Even though the value matrix of Table 1 only includes the value 
of the property of interest (porosity), the matrix may also contain the X, Y, Z location 
of the center point of the cell selected for the column and the area of the cell. The 
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matrix of Table 1 is derived by maximizing tlie porosity values, using a window that 
is three cell-layers thick, and filtering in a tight range (.236-.292 porosity); therefore, 
the numbers are fairly consistent and in the range of .75-.762 porosity (cumulative 
values of a three cell window). The data in Table 1 is fairly homogeneous; however, 
other datasets may have much more diversity. 



Table 1 





y1 


y2 


y3 


y4 


yS 


y6 


y7 


y8 


y9 


ylO 


y11 


y12 


y13 


y14 


y15 


y16 


y17 


y18 


x1 


0.492 


0.746 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


x2 


0.746 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


x3 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


x4 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


x5 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


x6 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


x7 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


x8 


0.761 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


x9 


0.761 


0.761 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


x10 


0.761 


0.761 


0.761 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


x11 


0.76 


0.761 


0.761 


0.761 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


x12 


0.76 


0.761 


0.761 


0.761 


0.761 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


x13 


0.76 


0.76 


0.761 


0.761 


0.761 


0.761 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


x14 


0.759 


0.759 


0.76 


0.76 


0.761 


0.761 


0.761 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 


0.762 



[039] Next in the method of one embodiment of the present invention, a 



summation matrix is generated based on the value matrix. In contrast to the value 
matrix, the summation matrix is used to determine the best cell location that's 
surrounded by other desirable cells (cells with desirable property values), not to 
determine the best single cell location. The summation matrix includes the 
accumulated values for each of the entries in the value matrix. Table 2 below 
illustrates an exemplary summation matrix, wherein the calculations described below 
have been applied to the value matrix of Table 1. 
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[040] In Table 2, the value in bold at location (XI 2, Y1 2) has been 
calculated to be the most desirable location of the processed three-dimensional 
model. 

Table 2 





y1 


y2 


y3 


y4 


y5 


y6 


y7 


y8 


y9 


ylO 


y11 


y12 


y13 


y14 


y15 


y16 


y17 


y18 


x1 


0.336 


0.37 


0.399 


0.422 


0.439 


0.448 


0.45 


0.448 


0.445 


0.441 


0.436 


0.429 


0.42 


0.407 


0.375 


0.349 


0.335 


0.279 


x2 


0.365 


0.389 


0.413 


0.433 


0.448 


0.456 


0.458 


0.456 


0.454 


0.45 


0.445 


0.438 


0.429 


0.414 


0.38 


0.358 


0.34 


0.278 


x3 


0.395 


0.414 


0.432 


0.447 


0.458 


0.465 


0.466 


0.465 


0.462 


0.459 


0.454 


0.448 


0.439 


0.422 


0.388 


0.368 


0.346 


0.276 


x4 


0.424 


0.438 


0.45 


0.461 


0.469 


0.473 


0.474 


0.474 


0.472 


0.468 


0.463 


0.457 


0.449 


0.43 


0.398 


0.379 


0.354 


0.275 


x5 


0 446 


0.457 


0.467 


0.473 


0.479 


0.482 


0.483 


0.483 


0.481 


0.478 


0.474 


0.467 


0.459 


0.438 


0.411 


0.392 


0.364 


0.275 


x6 


0.463 


0.472 


0.479 


0.484 


0.487 


0.489 


0.49 


0.49 


0.489 


0.486 


0.482 


0.476 


0.466 


0.445 


0.421 


0.402 


0.373 


0.276 


x7 


0.471 


0.477 


0.482 


0.485 


0.487 


0.49 


0.492 


0.494 


0.495 


0.496 


0.494 


0.49 


0.482 


0.456 


0.438 


0.419 


0.387 


0.28 


x8 


0.472 


0.479 


0.486 


0.492 


0.498 


0.503 


0.507 


0.511 


0.513 


0.514 


0.513 


0.51 


0.499 


0.475 


0.458 


0.436 


0.402 


0.283 


x9 


0.483 


0.491 


0.498 


0.504 


0.51 


0.515 


0.519 


0.523 


0.525 


0.526 


0.526 


0.524 


0.514 


0.492 


0.478 


0.452 


0.415 


0.297 


x10 


0.489 


0.496 


0.502 


0.509 


0.514 


0.519 


0.525 


0.531 


0.537 


0.542 


0.546 


0.548 


0.537 


0.52 


0.505 


0.478 


0.437 


0.31 


x11 


0.488 


0.497 


0.507 


0.517 


0.526 


0.534 


0.542 


0.55 


0.557 


0.563 


0.568 


0.571 


0.56 


0.545 


0.532 


0.503 


0.459 


0.316 


x12 


0.501 


0.511 


0.522 


0.531 


0.541 


0.549 


0.557 


0.565 


0.572 


0.578 


0.583 


0.586 


0.577 


0.565 


0.552 


0.523 


0.478 


0.336 


x13 


0 508 


0.517 


0.525 


0.533 


0.541 


0.548 


0.555 


0.561 


0.567 


0.573 


0.578 


0.582 


0.577 


0.572 


0.554 


0.536 


0.498 


0.358 


x14 


0.377 


0.383 


0.389 


0.395 


0.402 


0.409 


0.415 


0.422 


0.428 


0.434 


0.44 


0.446 


0.452 


0.45|0.448 


0.445 


0.426 


0.348 



[041] The summation matrix is populated by using the value matrix to 
calculate a distance-weighted sum of the values within a spacing radius of the center 
point of the cell selected in the value matrix. For each of the cells within the spacing 
radius selected by the user, the accumulated value (AccumValue) is determined 



using the relationships: 

[042] AccumValue = (CumWeightedValue/CuWeight) * CumWeightedArea, 

wherein 

[043] CumWelghtedValue = ^ value*weight, and 
[044] CumWeightedArea = Z area* weight, 

[045] Weight in the above calculations is derived using the relationship: 
weight = (SpacingRadius - DistanceFrom Cell)/SpacingRadius, wherein 
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DistanceFromCell is defined as the larger of the actual distance fronfi the cell for 
which the accumulated value is being calculated to the adjacent cell that is being 
taken into consideration or half the diagonal cell width (average distance to corner); 
and SpacingRadius is a user-defined value representing the resen/oir draining 
radius for each target. The SpacingRadius may be selected by the user based on a 
number of parameters (e.g., government regulations and/or financial 
considerations). 

[046] The calculations to arrive at the accumulated value are used to take 
into account the geometry of the various cells. If in the three-dimensional model, the 
cells are of relatively the same size then the above calculations may be 
unnecessary. However, in most models, the cells are of irregular geometry. For 
example, the center cells may be large cells, while the cells around the edges may 
be small cells. Therefore, if the cell values were summed or the cells were weighted 
based on distance, clusters of small cells may be Improperly favored over large 
cells. On the other hand, if the cells are weighted by cell area, large cells or cells 
that are barely within reach of large cells may be improperly favored over small cells. 

[047] In the accumulation calculations of one embodiment of the present 
invention, by taking into account the distance and the cell size, the irregularities of 
the cell geometry may be significantly reduced in the target selection process. 

[048] Once the summation matrix is complete, the summation matrix is 
analyzed to determine the most desirable value (0.586) in the matrix. (See Table 2 
above) Then, this value is selected as the first target location. Figure 7 provides an 
illustration of the location of the first target location on a three-dimensional model of 
the reservoir of Figures 2 and 3. The value in the matrix associated with the first 
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target and the value associate with other cells in the matrix that are within two times 
the user-defined radius are set to a value of zero for the selection of the next target 
location. 

[049] Table 3 below illustrates the result of the above described process, 
wherein the value of the selected target location and the values within two times the 
radius are set to zero. 



Table 3 





y1 


y2 


y3 


y4 


y5 


y6 


y7 


y8 


y9 


y10 


y11 


y12 


y13 


y14 


y15 


y16 


y17 


y18 


x1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.411 


0.399 


0.37 


0.342 


0.33 


0.282 


x2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.407 


0.375 


0.349 


0.335 


0.279 


x3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.38 


0.358 


0.34 


0.278 


x4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.368 


0.346 


0.276 


x5 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.379 


0.354 


0.275 


x6 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.364 


0.275 


x7 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.276 


x8 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.28 


x9 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.283 


x10 


0.483 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.297 


x11 


0.489 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


x12 


0.488 


0.497 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


x13 


0 501 


0.511 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


x14 


0.508 


0.517 


0.525 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



[050] Next from the remaining values in the summation matrix, the most 
desirable value (0.525) is selected as the second target location. (See Table 3 
above) Figure 8 provides an illustration of the location of the second target on a 
three-dimensional model of the reservoir displayed in Figures 2 and 3. In similar 
fashion to the process used to determine the next target after locating the first target, 
the value in the matrix associated with the second target and the value associate 
with other cells in the matrix that are within two times the user-defined radius is set 
to a value of zero for the selection of the next target. 
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[051] Table 4 below Illustrates the result of the above described process, 
wherein the value of the selected target and the values within two times the radius 
are set to zero. 



Table 4 







y2 


y3 


y4 


y5 


y6 


y7 


y8 


y9 


y10 


y11 


y12 


y13 


y14 


y15 


y16 


y17 


y18 


x1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.411 


0.399 


0.37 


0.342 


0.33 


0.282 


x2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.407 


0.375 


0.349 


0.335 


0.279 


x3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.38 


0.358 


0.34 


0.278 


x4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.368 


0.346 


0.276 


x5 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.379 


0.354 


0.275 


x6 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.364 


0.275 


x7 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.276 


x8 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.28 


x9 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.283 


x10 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.297 


x11 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


x12 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


x13 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


x14 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



[052] From the remaining values in Table 4, the most desirable value (.41 1 ) 



may be selected for the third target location. Then, the method repeats the above 
process of zeroing the value of the selected target location and the values within two 
times the radius of the selected target location to determine the next target location. 
The above process may continue until a user-defined maximum cutoff is reached or 
all possible targets have been located. The above examples use two times the 
radius; however, other values for the radius may be used. The above examples take 
into account one value of interest, but the method may be extended to use more 
than one value to determined the target locations. 

[053] The foregoing method may also be extended to take into account 
existing wells (e.g., exploratory wells). For example, the values in the summation 



17 



Attorney Docket No.: 2376.0016 

matrix within two times the radius of the existing well(s) may zeroed out before 
selecting the new targets. 

[054] Figure 9 provides an exemplary illustration of all the located targets 
(702, 802, and 902-952) on a three-dimensional model of the reservoir displayed in 
Figures 2 and 3. 

[055] Once all the possible targets have been located or the user-defined 
maximum has been reached, the target locations may be ranked, for example, from 
the most desirable (i.e., the first target location) to the least desirable (i.e., the last 
target location), and the user may display any percentage of the identified target 
locations. 

Methodology 

[056] Figure 10 is a flowchart illustrating the stages associated with the 
targeting method of one embodiment of the present invention. Method 1000 starts 
at (Stage 1001) and proceeds to Stage 1002. In Stage 1002, seismic data sets are 
used to build a three-dimensional model (Stage 1002) using an interpreter program. 
The program interprets the seismic data sets both for the top and bottom of the 
reservoir to develop a three-dimensional model from the interpreted structures which 
bound the reservoir as well as define the internal geometry. The three-dimensional 
model provides a representation of the characteristics of the reservoir underneath 
the surface. These characteristics, e.g., porosity, may be indicative of whether there 
are hydrocarbons in the reservoir and the fluidity of the hydrocarbons that may be 
found. 
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[057] Next, the three-dimensional model may be triangulated to create a 
three dimensional grid (X, Y, and Z axis) of cells (schematic model). (Stage 1004) 
Once the triangulation process is complete, for each X and Y coordinate, the 
schematic model is traversed in the Z direction to filter cell locations that do not meet 
a predetermined selection criteria (e.g., porosity > .2 and/or oil saturation > .3). 
(Stage 1006) 

[058] From the filtered three-dimensional model, the method of the 
preferred embodiment determines the most desirable set of contiguous cells for 
each X and Y location traversed in the Z direction (through the layers of the three- 
dimensional model) that meet one or more of the filtering criteria. The most 
desirable cells may be located using a windowing technique in which a sliding 
window composed, for example, of three cells or three layers may be used to find 
the location of the most desirable three continuous layers for each X and Y location. 
(Stage 1008) The window size may be determined by the number of layers in the 
three-dimensional model. 

[059] Then, the location of the center point or the midpoint of the most 
desirable contiguous cells for each X and Y location in the schematic model is 
identified, along with the associated value for the particular property of interest (e.g., 
porosity or oil saturation value), and entered into a two-dimensional matrix. (Stage 
1010) This process distills the schematic model down to a two dimensional value 
matrix representing the most desirable layers in the columns of the schematic 
model. Every column in the schematic model is represented by row and column 
indices (matrix entry) in the value matrix, and each matrix entry may contain the 
value of the property of interest (e.g., porosity or oil saturation). 
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[060] Next, a summation matrix is generated based on the value matrix. 
(Stage 1012) The summation matrix includes the accumulated values for each of 
the entries in the value matrix. The summation matrix is populated by using the 
value matrix to calculate a distance-weighted sum of the values within the spacing 
radius of the center point of the cell selected in the value matrix. (Stage 1012) For 
each of the cells within a spacing radius selected by the user, the accumulated value 
(AccumValue) is determined using the relationships: 

[061] AccumValue = (CumWeightedValue/CuWeight) * CumWeightedArea, 
wherein 

[062] CumWeightedValue = J] value*weight, and 

[063] CumWeightedArea = area* weight, 

[064] The "weight" in the above calculations may be derived using the 
relationship: weight = (SpacingRadius - DistanceFrom Cell)/SpacingRadius, wherein 
DistanceFromCell is defined as the larger of the actual distance from the cell for 
which the accumulated value is being calculated to the adjacent cell that is being 
taken into consideration or the half the diagonal cell width (average distance to 
corner); and SpacingRadius is a user-defined value representing the reservoir 
draining radius for each target. The SpacingRadius may be selected by the user 
based on a number of parameters (e.g., government regulations and/or financial 
considerations). 

[065] Once the summation matrix is complete, the summation matrix is 
analyzed to determine the most desirable value in the matrix. (Stage 1014) Then, 
the location associated with this value is selected as the first target location (Stage 
1016) and the value in the matrix associated with the first target and the value 
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associate with otiiers cell in the matrix that are within two times the user-defined 
radius are set to a value of zero for the selection of the next target location. (Stage 
1018) 

[066] Next, method 1000 proceeds to Stage 1020. In this stage, if a user- 
defined maximum cutoff for the number of targets has been reached or all the 
possible target have been located, method 1000 ends. (Stage 1020 and 1022). 
However, if the above is not true, method 1000 returns to Stage 1016 and selects, 
from the remaining values in the summation matrix, the most desirable value as the 
next target location. Then, in similar fashion to the process used to determine the 
previous target, the value in the matrix associated with the current target and the 
value associate with other cells in the matrix that are within two times the user- 
defined radius are set to a value of zero. 

[067] Again, method 1000 returns to Stage 1020 and determines if a user- 
defined maximum cutoff for the number of targets has been reached. If so, method 
1000 ends. Otherwise, the method returns to Stage 1016 and 1018 for the selection 
of the next target. The loop created by Stages 1016-1020 may be executed until the 
cutoff is reached or all targets have been located; then method 1000 ends. (Stage 
1022) 

[068] Other embodiments of the invention will be apparent to those skilled 
in the art from consideration of the specification and practice of the invention 
disclosed herein. It is intended that the specification and examples be considered 
as exemplary only, with a true scope and spirit of the invention being indicated by 
the following claims. 
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